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Direct numerical simulations of electron dynamics in externally driven electrostatic waves have 
been carried out using a relativistic two-fluid one-dimensional Vlasov-Poisson code. When the driver 
wave has sufficiently large amplitude, ion density holes (cavitons) form. The interaction between 
these cavitons and other incoming Langmuir waves gives rise to substantiai local acceleration of 
groups of electrons, and fine jet-like structures arise in electron phase space. We show that these 
jets are caused by wave-breaking when finite amplitude Langmuir waves experience the ion density 
gradient at the leading edge of the holes, and are not caused by caviton burn-out. An analytical two- 
fluid model gives the critical density gradient and caviton depth for which this process can occur. In 
particular, the density gradient critically affects the rate at which a Langmuir wave, moving into the 
caviton, undergoes Landau damping. This treatment also enables us to derive analytical estimates 
for the maximum energy of accelerated electrons, and for the energy spectrum along a phase-space 
jet. These are confirmed by direct numerical simulations. 

PACS numbers: 52.35.Mw, 52.38.Kd, 56.65.Ff 



I. INTRODUCTION 



Outline 



This paper describes a novel mechanism for the non- 
resonant acceleration of electrons within a plasma. We 
demonstrate how coupling of strong Langmuir turbulence 
to finite amplitude, coherent Langmuir waves allows the 
potential energy of the Langmuir waves to be imparted to 
electrons as the wave breaks. This results in the accelera- 
tion of a small population of fast electrons from the back- 
ground population at the point where wave-breaking oc- 
curs. Acceleration occurs spontaneously without requir- 
ing an intense laser pulse or particle injection, which are 
key features of electron acceleration schemes QlJlJlJlI 
that rely on wave-particle resonance, such as laser wake- 
field acceleration (LWFA), self- modulated LWFA, plasma 
beat wave accelerators and plasma wakefield accelerators. 
The electrons accelerated by our new non-resonant mech- 
anism form a collimatcd phase space jet in the electron 
distribution function f e . Such populations of energetic 
electrons are of some concern in inertial confinement fu- 
sion (1CF) experiments, because they may contribute di- 
rectly to capsule pre-heat or form the seed populations 
for further acceleration. For example, the hot electrons 
generated by this mechanism could undergo wakefield ac- 
celeration m the presence of a laser field, creating 
even more ener geti c particles capable of pre-heating the 
1CF capsule Hlli3. 

By considering the relativistic one dimensional Vlasov- 



Poisson system and exposing this system to a large am- 
plitude electrostatic standing wave (in a similar manner 
to Ref . [TH| ) . it is possible to drive the formation of cavi- 
tons. These cavitons are, in essence, local regions of ion 
density depletion accompanied by a local concentration of 
electric field amplitude. Their formation is driven by the 
ponderomotive force and they are seeded by small fluc- 
tuations in the electron density. In our simulations these 
initial fluctuations are provided by the electrostatic driv- 
ing field. However, once the process of cavitation has 
begun, the external field can be removed without sig- 
nificantly affecting further caviton development. In this 
paper we show that at a particular point in the caviton 
development, fine jet-like structures form in the electron 
phase space, and we derive the necessary conditions for 
jet formation depending on the characteristics of the cavi- 
ton. These phase space jets represent a portion of the 
electron population accelerated from a small region near 
the edge of a caviton. The process by which these the jets 
emerge, together with their characteristics, suggests that 
they are the result of breaking Langmuir waves. This 
hypothesis is supported by further simulations, and by a 
simple model for the process based on a fluid treatment 
including Landau damping. This theor y is an extension 
of the work of Akhiezer and Polovin [12| on breaking 
Langmuir waves in a uniform medium. The energy dis- 
tribution and the maximum energy within a jet can be 
calculated from conservation of energy, provided we as- 
sume that all of the potential energy of the wave is im- 
parted to the electrons at the breaking wavefront. 
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B. The Relativistic Vlasov-Poisson System 

The model used in this paper is a one dimensional rel- 
ativistic Vlasov-Poisson system of electrons and protons 
with a mass ratio mi/m e = M r , and no magnetic field. 
This fully nonlinear self consistent system is governed by 
the Vlasov equation for the electron distribution function 
fe(x,p,t) 

2Il + _ eE— = (1) 

dt m e j dx dp 

the Vlasov equation for the ion distribution function 

fi(x,P,t) 

d -£ + JL f+eEf = 0, (2) 
Ot rriij ox op 

and Poisson's equation for the electric field 

i— £(/'■* -/**) <3) 

where x is the spatial co-ordinate, p = ^yin e v is the 
momentum co-ordinate and 7 = (1 + p 2 /m^c 2 ) 1 ' 2 is 
the Lorentz factor. Taking z to represent the nor- 
malised form of the variable z, the following dimension- 
less normalisation, appropriate to a relativistic system, is 
adopted throughout: x — (c/ui pe )x, t — (l/uj pe )t, v = cv, 
E = {u> pe cm e /e)E, p — m e cp. It follows that frequen- 
cies are normalised to the plasma frequency, uj = Lo pe LO, 
wavenumbers to the ratio of the plasma frequency to the 
speed of light, k — uj pe k/c, and temperatures are nor- 
malised such that T e .i = (kB/m ei iC 2 )T ei i, where fee is 
Boltzmann's constant. All simulations are carried out us- 
ing a mass ratio M r — 100. This value is sufficiently large 
to allow the development of phenomena on two disparate 
time-scales (since electron and ion plasma frequencies are 
an order of magnitude apart) , without the increased run- 
time of the real ratio. 



C. Numerical Approach 



to the high quiver velocity regime v 2 /v 2 -~ e > 1, where 
v q = eE /m e oj Q . Hence E Q > m e u> VT e /e, or in nor- 
malised units 

E > VTe (4) 

This driving field, similar to the perturbation used in ear- 
lier work on the modulational instability , is necessary 
to drive the formation of cavitons and is only needed dur- 
ing the early stages of the simulation. In our simulations 
it is applied from f=0tof=10, which is sufficient to 
seed the formation of cavitons. 

Extending the presence of the driving field after t = 10 
does not significantly effect the evolution of cavitions or 
the amplitude of Langmuir waves in the system. Hence 
the energy of accelerated electrons remains < 5MeV. 
However, Langmuir waves formed as a result of the the 
continued presence of E4 result in more numerous and 
frequent phase space jets. 

II. RESULTS 

A. Initial Conditions 

The system is initialised with a Maxwellian distribu- 
tion of both electrons and ions at equal temperatures 
with an electron thermal velocity 

v 2 e = c 2 /10 3 (5) 

and an ion thermal velocity 

v 2 Tl = c 2 /10 3 M r , (6) 

equivalent to a temperature ~ 0.5keV. In normalised 
units this gives vx e = 10~ 3 / 2 , hence Eq.QJ implies that 
we require Eq > 10~ 3 / 2 . The simulation box is of length 
L x = 207r in normalised units, with periodic boundary 
conditions in space. In summary, the dimensionless ini- 
tial conditions are f e = 0.001, % = 0.001, L x = 20tt, 
Eq = 0.5, cJo = 1, k — 2ir/L x , and the initial maximum 
of f e , fr ax = 9.2925. 



The Vlasov-Poisson system is solved using the code 
detailed in Ref.^3]. This is a split Eulerian scheme in 
which the distribution functions (/ e , ft) are calculated on 
a fixed Eulerian grid, and the solver is split into separate 
spatial and velocity space updates Q. These updates 
are one dimensional, constant velocity advections car- 
ried out using the piecewise parabolic method [l5|. The 
original code has been extended to solve the fully rel- 
ativistic Vlasov-Poisson system Eqs.Q to |J3J. A large 
amplitude external driving field Ed = Eo sin(fci) sm(ujat) 
is added to the self consistent electric field found from 
Poisson's equation. The parameters of the driving field 
are chosen to ensure that the system is driven at res- 
onance (loq = Lo pe , or ljq = 1 in normalised units), 
and that the intensity of the perturbations corresponds 



B. Caviton Formation and Jet Creation 

The development of cavitons can be seen in both the 
isometric renderings of the electron distribution function 
(Figs and [2J and plots of the density and electric field 
amplitude (Fig|3J) . Cavitons form in response to the pon- 
deromotive force exerted on the electrons by the driving 
field Ed- Hence they are centred around regions where 
the amplitude of the driving field is strongest, leading to 
symmetry in n e ^ and \E\ about the centre of the system 
x = lOir. By t = 20 two prominent holes, driven by the 
ponderomotive force, have formed in the electron density 
centred at x « 25 and 37. As the system continues to 
evolve, the ions are slowly evacuated to form two broad 
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FIG. 1: Isometric view of the electron distribution function f e (left) and contour plot o/log(/ e > 10~ 6 ) (right) at time t = 20. 
Here the relativistic Vlasov-Poisson system is driven from t — to 10 at u) — ui pe . The evacuation of electrons from localised 
regions of the simulation domain is seen in the electron distribution function; this is the early stage of cavitation. Caviton 
formation is due to the modulational instability initially driven by the external field Ed ■ 




FIG. 2: Isometric view of the electron distribution function f e (left) and contour plot o/log(/ e > 10 -6 ) (right) at time t = 40. 
Here the relativistic Vlasov-Poisson system was driven from t — to 10 at ui = ui pe . Fine jet-like structures in electron phase 
space are visible: jets marked 'a' have recently formed on the inside edges of the two cavitons, whereas the jets marked 'b' have 
formed in the same area at an earlier time and have since been advected through the system. The jets marked 'c ' are at an 
intermediate stage, having formed on the outside edges of the cavitons. These jets are the result of Langmuir wave-breaking at 
the edges of the evolving density holes. 



density holes. These holes in both ion and electron den- 
sities are accompanied by a local increase in the electric 
field amplitude, characteristic of caviton formation, see 
FigEl 

At later times, populations of accelerated electrons are 
clearly visible in the contour plots of the electron dis- 
tribution function. These are the electron phase space 
jets. At t = 20, ten plasma periods after the driving field 
has been removed, there is no evidence of jet formation, 
see Fig^ However, between i — 20 and t — 40 a series 
of phase space jets form on both edges of the deepening 
cavitons. Studying the evolution of the electron distri- 
bution function during this period allows one to identify 
phase space jets at different points in their evolution. In 



Fig0we highlight jets at three separate stages of evolu- 
tion. There are two jets (a) forming on the inside edges of 
the caviton as well as two old jets (b) which have been ad- 
vected across the system, effectively crossing the cavitons 
on whose edges they formed. Finally, there are two inter- 
mediate jets (c), on the outer edges of the cavitons. The 
most energetic electrons within these jets have energies s=s 
5MeV. The appearance of these jets does not affect the 
development of the cavitons, which continue to deepen 
after the appearance of the phase space jets. In addition, 
the simulation is seeded with randomly distributed tracer 
particles. These do not contribute to the numerical solu- 
tion of the Vlasov Poisson system and are simply moved 
in response to the total electric field (self consistent elec- 



4 




■ 



i.o 



'-5 



5 



A f\ 



Kmj 

■ 

- 



C. Jet Emergence 

The process of jet emergence is best explained by ob- 
serving the evolution of the electron distribution func- 
tion, focusing on a region where a jet develops, during 
its early stages. Figure 0] shows a reduced section of the 
electron phase space, the region 50 < x < 60 where the 
rightmost jet, labelled 'b' in FigJ^J first appears. This re- 
gion encompasses the right hand side of one of the deep- 
ening cavitons seen in Fig|3| The jet forms at the outer, 
right hand edge of this caviton at x « 60. The sequence 
of contour plots in FigQ] shows that the jet then extends 
out from the main electron distribution, at \p\ « 2, to 
momenta of \p\ ~ 8 in approximately one plasma period. 



FIG. 3: Plots of the electron (solid line) and ion (dashed 
line) density at t = 20 (top left) and t — 40 (top right), 
together with the electric field amplitude at t = 20 (bottom 
left) and t = 40 (bottom right) for a relativistic Vlasov- 
Poisson system that was driven from t = to 10 at to — 
ujpe . The formation of ion density holes can be clearly seen, 
together with localised growth in the electric field amplitude 
characteristic of caviton formation. Caviton formation is due 
to the modulational instability driven initially by the external 
field Ed- 
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FIG. 4: Contour plots o/log(/ e > 10~ 6 ) in the region 50 < 
x < 60 at times t — 32 (top left), 34 (top right), 36 (bottom 
left) and 38 (bottom right). These plots show a phase space 
jet developing at the outside edge of the caviton (see Fig\^j 
at x » 60. Electrons are accelerated from a compact region 
of the background distribution to form the jet, which extends 
and advects across the caviton. 



trie field plus driving field) . Following the motion of these 
particles indicates that the electrons which form the jets 
are not accelerated from within the caviton - this effec- 
tively rules out caviton burn-out as an explanation for 
the origin of the phase space jets. 



III. PHYSICS OF ELECTRON JETS 
A. Wave Breaking 

The key observed features of the phase space jets are: 

1. Jets do not appear at early times in the simulation, 
they first require some degree of caviton evolution. 

2. Jets are not directly related to the external driver, 
they appear after the removal of the external driv- 
ing field. 

3. Jets are not associated with caviton burn-out pro- 
cesses, since cavitons persist long after the appear- 
ance of jets. 

4. Jets originate at the caviton edge. Electrons are 
accelerated from the main distribution to form the 
phase space jet which is then advected through the 
system. Their constituent electrons pass over the 
caviton on whose edge they formed, escaping the 
influence of the caviton completely. The direction 
of this advection indicates that phase space jets are 
the result of processes originating outside the cavi- 
ton. 

5. Electrons are accelerated up to energies of 5MeV 
from an initially Maxwellian population with tem- 
perature w 0.5keV. 

From extensive numerical simulations of the system, it is 
clear that the breaking of Langmuir waves on the density 
gradients at the edges of the cavitons is responsible for 
the creation of the phase space jets. Figure [S] provides a 
schematic illustration of the physical process, which we 
explore in the rest of this section. The Langmuir wave 
first approaches the density hole. As it moves into the 
region of lower density, the phase speed at the front of 
the wave falls. If this proceeds rapidly enough to over- 
come the effect of Landau damping (which acts to damp 
the incoming wave energy, and thereby prevent it from 
breaking), then the wave will break. This creates a strong 
electric field localised at the wave crest, which accelerates 
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FIG. 5: Schematic representation of electron acceleration by 
Langmuir wave-breaking, a: Density hole with a few sample 
electrons highlighted, b: The crest of a Langmuir wave moves 
towards the density hole, c: As the wavefront moves down the 
gradient at the edge of the density hole, its phase velocity falls 
and the wavefront steepens, d: For a density hole of suffi- 
cient depth and steepness, the incoming wave breaks, forming 
a strong localised electric field which accelerates electrons from 
the inside edge of the caviton to create jets in electron phase 
space. 



electrons in the vicinity away from the background popu- 
lation to form phase space jets. In the next three sections 
we outline an analytical model of the competing processes 
of wave-breaking and Landau damping, together with a 
derivation of the electron energy distribution within the 
jet. 



B. Wave Breaking Amplitude 

The maximum amplitude Eb r of a Langmuir wave 
which can be sustained in a plasma before it breaks is 
given by the Akhiezer and Polovin constraint For a 
non-relativistic phase velocity v p , this is given by 



eE br /m e ujp e v p = 1 



(7) 



where the breaking amplitude Ebr is the maximum elec- 
tric field amplitude which can be sustained by a Lang- 
muir wave before it breaks. In normalised units Eq.Q) is 
equivalent to Ey. — v p ; the breaking amplitude depends 
on the phase velocity of the wave, which is determined 
by the Langmuir dispersion relation 

uj 2 = Lu 2 pe + 3v Te 2 k 2 (8) 

Consider a Langmuir wave moving through a region of 
decreasing density 

ft = n(x) (9) 
where n falls from an initial value no, corresponding to 

1/2 

a plasma frequency of ojp e = (noe 2 /m e eo) ■ As a first 
approximation, we may assume that the amplitude and 
frequency of the Langmuir wave remain fixed at El (El 
in normalised units) and ui « uP pe (Q sa 1 in normalised 
units) respectively. Taking the dispersion relation Eq. JSJ 
in normalised units 



Q 2 



n(x) + 3T e k 2 



(10) 



with w « 1 gives an estimate of the local wavenumber of 
the Langmuir wave as a function of n(x): 



k = 



((1 - h(x)) /3T e 



~\l/2 



(11) 



This corresponds to a phase velocity v p — & jk w 1/fc 
given by 



v p = (3T e /(l - n(x)) 



1/2 



= E, 



br 



(12) 



the breaking amplitude, by Eq.Q. For the case of a 
linear density ramp n(x) = 1 — i}x where < x < l/n, 
Eas. lfTT|l and ljT2|) imply 



- 2\ V2 

k = [r\xj 3T e 



and 



v p = \ 3f e /r]x\ 



1/2 



Ebr 



(13) 



(14) 



As the wavefront moves down the density ramp, the 
phase velocity falls, reducing the maximum wave ampli- 
tude which can be sustained. If the breaking threshold 
falls sufficiently, it will be satisfied by the incoming wave, 
which will then break. The wave breaking condition 



Ebr < El 
together with Eq. l|13fl gives 

hT e / V x) <E L 



(15) 



(16) 



Thus far, this treatment has not taken into account the 
effect of Landau damping of the Langmuir wave as it 
moves down the density ramp, which could reduce the 
field amplitude at a rate which ensures the breaking con- 
dition is never achieved. If the initial amplitude of the 
Langmuir wave is Elo, then the damped amplitude at 
time t is given by 



Er 



E L oe~ r V 



(17) 



Here the local damping decrement jl determines 

r(t) 



r (t) = / lL dt 



1L 



8/ fc 3 A^ 



exp 



2 2k 2 \ 2 D 



(18) 



(19) 



Using the estimate for k given in Eq. (|13fl and solving the 
integral in Eq. l|18|) for x following the wave, given that 



dx 
~dt 



dui 
~dk 



\3f e r)xJ 



1/2 



(20) 



yields 

^ = ^ Loexp (^k©" cxp H^)) (21) 

Combining Eas. l21() and (|15fl gives the condition for 
wave-breaking, and hence phase-space jet formation 

^° exp (^k© %xp H^)) 

-1^1 > (22) 
\r)X J 

for a Langmuir wave with initial electric field amplitude 
Elq moving down a linear density ramp with gradient 
—Tj. In the steep gradient limit (i.e. 77 — ► 00), this is 
equivalent to Ea. (|16l) . To summarise, as the Langmuir 
wave moves down the density ramp it experiences Lan- 
dau damping which reduces the wave amplitude from its 
initial value. The wave-breaking condition is then only 
satisfied if the decline in density reduces Eb r to the value 
where it is equal to the damped Langmuir wave ampli- 
tude El given by Ea. l21|) . Thus for a density ramp with 
edge gradient 77 and depth A, the wave will break, thereby 
accelerating electrons to form phase space jets, if Eq.J32J) 
is satisfied for some x, < x < A/77. This relation has 
been derived from a simple fluid treatment, modified to 
take into account the effect of Landau damping, which 
depends on caviton gradient 77 and depth A. It is sup- 
ported by numerical simulations using fixed background 
density profiles containing cavitons with a chosen depth 
and edge gradient, which are presented below. 
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FIG. 6: Top: Schematic of the initial density profile used for 
fixed ion simulations. The simulation boundaries are periodic 
with density holes defined by a 'tanh' function. A traveling 
Langmuir wave of fixed wavelength and amplitude is added 
in the form of a sinusoidal electron density and velocity per- 
turbation at the centre of the system. Bottom: Isometric 
rendering of the electron distribution function from the fixed 
ion simulations at t = 0. 



C. Fixed Ion Simulations 

In order to test the model outlined above and demon- 
strate that the phase space jets are a direct result of 
breaking Langmuir waves, a further series of numerical 
simulations was conducted. These simulate the interac- 
tion of a Langmuir wave with a system having immobile 
ions whose initial density profile is fixed to provide den- 
sity holes of a prescribed size and shape, as shown in 
Fig|S] A Langmuir wave was created by perturbing the 
initial electron distribution function to create a traveling 
wave of a given amplitude and wavelength. This initial- 
isation relies on the linear dispersion relation, and so is 
only valid for small amplitudes. As the wave propagates 
across the simulation domain, it encounters the density 
hole created by the fixed ion background density profile. 
By varying the parameters of the density hole and the 
Langmuir wave, it is possible to change the nature of the 
phase space jets and confirm the functional dependence 
on caviton depth and gradient predicted by Eq. 1221) . Fig- 
ure0shows an area of (A, 77) parameter space partitioned 
by Ea. (|22() . for a Langmuir wave of amplitude Elq = 0.12 
and wavelength Al = 2tt. It also displays the results of 



simulations of the interaction of this Langmuir wave for 
a range of (77, A), classified by the presence or absence of 
phase space jets. Our criterion for jet formation is the 
presence of a second maximum, separated from the main 
distribution, in the region of the density ramp. This re- 
quires that f e pass through a critical value f c four times 
for some point in the density ramp. The value f c is chosen 
to be 1% of the initial maximum f™ ax of f e . The sim- 
ulation results support the analytical breaking condition 
Ed. lf2~2|l derived above. Figure |S| shows such a simulation 
at t = 56: a series of phase space jets have formed as 
a Langmuir wave (Xl = 27r, Sn e = 0.12) encountered a 
density hole of depth A = 0.6 and gradient 77 = 0.5. 



D. Electron Energy Distribution within Jets 

The maximum energy achieved by electrons in a phase 
space jet, and the distribution of electrons within the 
jet, can be calculated from the assumption that the total 
energy carried by the sum of jet electrons comes from 
the potential of the breaking Langmuir wave. Here we 
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FIG. 7: Jet formation and non-formation, plotted with re- 
spect to caviton parameter space (A, rf), where A represents 
the caviton depth and n the edge gradient. The space is parti- 
tioned into regions of jet formation and non-formation by the 
breaking condition Ea. and superimposed are the results 
from simulations using a fixed ion background to produce the 
desired caviton structure. The Langmuir wave considered here 
has amplitude Elo — 0.12 and wavelength Al = 2ty. 



consider the fully relativistic case and show that conser- 
vation of energy enables us to find the electron energy 
distribution along the jet. Following the same approach 
as Ref.0], we assume the energy gain of electrons ac- 
celerated by the breaking wave is a positive, continuous, 
single-valued function of their initial position (i.e. their 
proximity to the breaking wavefront). We then construct 
an expression for the electron density as a function of 
energy gain Ae. However, we must account for the pos- 
sibility that two electrons from different spatial positions 
could achieve the same energy gain. We therefore parti- 
tion the domain into regions where the number density 
of accelerated electrons is a single valued function of Ae, 
as shown schematically in Fig[§] which we discuss below. 
The total number of accelerated electrons N Jet is given 
by integrating the electron density over the spatial region 
So from which electrons are accelerated: 



N Jet _ 



So 



n(xo)dxo 



(23) 



where So represents the initial position of the electron. 
This can be rewritten in terms of a new variable Ae(xo), 
the energy gain as a function of initial position, to give 



N 



Jet 



E 



rij (x (Ae)) 



dAe 



dxo 



dAe 



(24) 



Here the region of integration So is partitioned into sub- 
domains {Qj}jeN- within each sub-domain Qj, the func- 
tion Ae(io) is single valued and has no turning points, 
as illustrated in Fig|21 The boundaries between the sub- 
domains are defined by the location of the turning points 
of Ae(aio)- Denoting the energy spectrum within the jet 
as a function of energy gain by n Jet (Ae), we can write the 
total number of electrons in the jet N Jet as an integral 



of n Jet over Ae 



N Jet = I n Jet( Ae } dAe 

So 



Equating this to Ea. (|24() gives 

et (Ae)=5>i(a:o (Ae)) 



, Jet 



dAe 



dxo 



(25) 



(26) 



The function xo(Ae) comes from the potential function 
of the Langmuir wave. Since the electric field of the wave 
near to breaking will be of nonlinear sawtooth form, its 
potential <f>L can be represented 01 by a parabolic func- 
tion with amplitude 4>o 



4xg 



(27) 



over the range of one wavelength Al. From the as- 
sumption that electrons accelerate using the local po- 
tential energy of the Langmuir wave, it follows that 
Ae(io) = 4>l{xo), with the maximum possible energy 
gain Ae max = 4>o- 



Ae = Ae r , 



Inverting this expression yields 




(28) 



xo = Xl ( 1 + (1 + Ae/Ae max ) 1/2 ) /2 (29) 



while differentiating with respect to Xq gives 



dAe 



dxo 



= Ae 



SXo 



(30) 



Substituting Eas.l|2*^)l and (JSOJ into Ea.l|2l>|). for a con- 
stant background density profile, gives the energy spec- 
trum for a single phase space jet: 



(Ae) 



, Jet 



(0) 



1 



(1 + Ae/ Ae ma x) 



1/2 



(31) 



Here we have normalised the equation using n Jet {Ae = 
0) = no\T,/^Ae maxi for a local electron density of no- It 
follows that the ratio of highest energy to lowest energy 
electrons within the jet is given by 



Jet 



n— (Ae max )/n Jet {Ae = 0) = 1/V2 



(32) 



The analytically derived energy spectrum Eq. (|31|l , shown 
in FigE3 (left), is broadly consistent with the electron 
distribution function that is obtained from the simula- 
tions discussed earlier. Figure ^| (right) shows a plot of 
f e , normalised to its value f e o at the base of the phase 
space jet, along the ridge of the phase space jet shown in 
Fig|l' 
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FIG. 8: Isometric view of the electron distribution function f e (left) and contour plot o/log(/ e > 10 6 ) (right) at time t — 56 
from a simulation of the interaction of a Langmuir wave with an ion density hole ( created and maintained via the prescribed 
immobile ion distribution of Fig\^ top). A Langmuir wave (Xl = "2n, 8n e — 0.12) has encountered the density hole of depth 
A = 0.6 and gradient r\ = 0.5. The incoming wave breaks as it moves into the density hole, producing a strong localised electric 
field which accelerates electrons from the bulk of the distribution to form a jet. 




x o a&e/fn = n x 



FIG. 9: Schematic representation of the decomposition of the 
spatial domain So into sub-domains {Qj}jsn at equation \2$ . 
The domain So represents the region of space from which elec- 
trons are accelerated to form a phase space jet. In order to 
find the energy spectrum, n(Ae) for a phase space jet, we re- 
quire Ae(x). However, we must account for the possibility 
that Ae could be multi-valued with respect to x, for example 
Ae(xa) = Ae(xb) = Ae a> 6. Dividing So into {Qj}je!<t, where 
the sub-domain boundaries are defined by the condition that 
dAe/dx = 0, avoids the problem of Ae(x) being multi-valued. 
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IV. CONCLUSIONS 

Shortly after exposing a fully relativistic Vlasov- 
Poisson system to a strong external driving field at res- 
onance ujq = uj pe , we observe the formation of fine struc- 
tures in the electron distribution function, corresponding 
to the acceleration of small populations of electrons to 
high energies (~ 5MeV). These phase space jets result 
from the breaking of Langmuir waves (initially excited 
by the driving field) in the density gradients of cavitons 
formed via the ponderomotive force exerted by the exter- 
nal driving field. Restricting the presence of the driving 
field to t < 10 allows direct acceleration by the driving 
filed to be ruled out as a mechanism. We are able to 
explain the origin of the small populations of electrons 
which form the jets; derive a necessary condition for the 
formation of jets starting from a basic fluid treatment, 
which is supported by further numerical simulations; and 
estimate the energy distribution of electrons within the 
phase space jet. This process may arise whenever Lang- 
muir waves move through a strong density gradient, and 
is not limited to one dimension or to caviton formation. 
It may therefore require consideration in laser-plasma in- 
teraction contexts spanning incrtial confinement fusion 
and particle acceleration. 



FIG. 10: Similarity of predicted and measured forms of the 
energy distribution of electrons within a jet. Left: Analytical 
result Eq. 13 IV . Right: The value of the distribution function 
f e , normalised to its value f e o at the base of the phase space 
jet. The energy gain of electrons in the phase space jet is 
normalised to the maximum observed energy ~ 5MeV. 
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